Multi-Higgs U(l) Lattice Gauge Theory in Three Dimensions 
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We study the three-dimensional compact U(l) lattice gauge theory with TV Higgs fields numeri- 
cally. This model is relevant to multi-component superconductors, antiferromagnetic spin systems 
in easy plane, inflational cosmology, etc. For TV — 2, the system has a second-order phase transition 
line 61(02) in the ci (gauge coupling)— ci (Higgs coupling) plane, which separates the confinement 
phase and the Higgs phase. For TV = 3, the critical line is separated into two parts; one for C2 < 2.25 
with first-order transitions, and the other for C2 > 2.25 with second-order transitions. 
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There are many interesting physical systems involving 
multi-component (TV-component) matter fields. Some- 
times they are associated with exact or approximate sym- 
metries like "flavor" symmetry. In some cases, the large- 
TV analysis is applicable and it gives us useful informa- 
tion. But the properties of the large- TV systems may dif- 
fer from those at medium values of TV that one actually 
wants to know. Study of the TV-dependence of various 
systems is certainly interesting but not examined well. 

Among these "flavor" physics, the effect of matter 
fields upon gauge dynamics is of quite general inter- 
est in quantum chromodynamics, strongly correlated 
electron systems, quantum spins, etc.0] In this let- 
ter, we shall study the three-dimensional (3D) U(l) 
gauge theory with multi-component Higgs fields <f> a (x) = 
|</> a (aOI exp(iv? a (x)) (a = 1, • • ■ , TV). This model is of gen- 
eral interest, and knowledge of its phase structure, order 
of its phase transitions, etc. may be useful to get better 
understanding of various physical systems. These sys- 
tems include the following: 

TV -component superconductor: BabaevQ argued that 
under a high pressure and at low temperatures hydro- 
gen gas may become a liquid and exhibits a transition 
from a superfluid to a superconductor. There are two 
order parameters; <p e for electron pairs and <f) p for pro- 
ton pairs. They may be treated as two complex Higgs 
fields (TV = 2). In the superconducting phase, both cf> e 
and ipp develop an off-diagonal long-range order, while 
in the superfluid phase, only the neutral order survives; 
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p-wave superconductivity of cold Fermi gas: Each 
fermion pair in a p-wave superconductor has angular mo- 
mentum J = 1 and the order parameter has three com- 
ponents, J z = —1,0,1. They are regarded as three Higgs 
fields (TV = 3). As the strength of attractive force be- 
tween fermions is increased, a crossover from a super- 
conductor of the BCS type to the type of Bose-Einstein 
condensation is expected to take place 4]. 

Phase transition of 2D antiferromagnetic (AF) spin 
models: In the s = 1/2 AF spin models, a phase tran- 



sition occurs from the Neel state to the valence-bond 
solid state as parameters are varied. Senthil et al.[|[ 
argued that the effective theory describing this transi- 
tion take a form of U(l) gauge theory of spinon (CP 1 ) 
field z a (x) [\zi\ 2 + \z2\ 2 = 1). In the easy-plane limit 
(S z = 0), \zi\ 2 = \z 2 \ 2 = 1/2 and so they are expressed 
by two Higgs fields as z a — exp(itp a )/\/2 (TV = 2)0. 

Effects of doped fermionic holes (holons) to this AF 
spins are also studied extensively. The effective theory 
obtained by integrating out holon variables may be a 
U(l) gauge theory with TV = 2 Higgs fields (with non- 
local gauge interactions). Kaul et al.pj predicts that 
such a system exhibits a second-order transition, while 
numerical simulations of Kuklov et al. Q exhibit a weak 
first-order transition. This point should be clarified in 
future study. 

Inflational cosmology: In the inflational cosmology [9|, 
a set of Higgs fields is introduced to describe a phase 
transition and inflation in early universe. Plural Higgs 
fields are necessary in a realistic modelficj. 

The following simple consideration "predicts" the 
phase structure of the system. Among TV phases ip a {x) 
of the Higgs fields, the sum ip + = VJ a <Pa couples to the 
gauge field and describes charged excitations, whereas 
the remaining TV — 1 independent linear combinations 
(f>i(i = 1, • • • , TV — 1) describe neutral excitations. The 
latter TV — 1 modes may be regarded as a set of N — 1 
XY spin models. As the TV = 1 compact U(l) Higgs 
model stays always in the confinement phase [llj. we ex- 
pect N — 1 second-order transitions of the type of the 
XY model. 

Smiseth et al.[12j studied the noncompact U(l) Higgs 
models. A duality transformation maps the charged sec- 
tor into the inverted XY spin model. Thus they pre- 
dicted that the system exhibits a single inverted XY 
transition and N — 1 XY transitions. Their numerical 
study confirmed this prediction for TV = 2. 

For TV = 2, Kragset et al.[l3[ studied the effect 
of Berry's phase term in the TV = 2 compact Higgs 
model. They reported that Berry's phase term sup- 
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FIG. 2: Instanton density p for N 
function of ci. 
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FIG. 1: (a) System-size dependence of specific heat C for 
N = 2 at C2 = 0.4. (b) Scaling function 7j(x) for FigQTa). 



presses monopoles (instantons) and changes the second- 
order phase transitions to first-order ones. 

In this letter, we shall study the multi-Higgs models 
by Monte Carlo simulations. We consider the simplest 
form, i.e., the 3D compact lattice gauge theory without 
Berry's phase; the Higgs fields <f> xa are treated in the 
London limit, \<fi xa \ = 1. The action S consists of the 
Higgs coupling with its coefficient c\ a {a = 1, . . . , N) and 
the plaquette term with its coefficient C2, 

1 N 

S = \Cla(t>l + ^ a Uxn^xa + H.C.) 

+y Y.^ u l+u, iJ P x+ll , v u xll + B.. c ), (l) 

where U xfl [= exp(i9 xll )} is the compact U(l) gauge field, 
fi, v{= 1, 2, 3) are direction indices (we use them also as 
the unit vectors). 

We first study the N = 2 case with symmetric cou- 
plings Cn = c 12 = C\. We measured the internal energy 
E = -(S)/L 3 and the specific heat C = ((S - (S)) 2 ) / L 3 
in order to obtain the phase diagram and determine the 
order of phase transitions, where L 3 is the size of the 
cubic lattice with the periodic boundary condition. 

In Fig[lja), we show C at C2 = 0.4 as a function 
of c\. The peak of C develops as the system size 
is increased. The results indicate that a second-order 
phase transition occurs at c\ ~ 0.91. By applying the 
finite-size-scaling (FSS) hypothesis to C in the form of 
C(d,L) = L'^L^e), where e = (a - c loo )/c loo and 
cioo is the critical coupling at L — * oo, we obtained 
v = 0.67, cr = 0.16, and ci^ = 0.909. In FigQIb) we 
plot r](x), which supports the FSS. 

The above results for N — 2 are consistent with the 
"prediction" given above. The sum (p x+ = tp x i +ip X 2 cou- 
ples with the compact gauge field and generates no phase 
transition ll|, while the difference tp x - = ip x \ — Lp x i be- 
haves like the angle variable in the 3D XY model. The 
3D XY model has a second-order phase transition with 
the critical exponent v = 0.666...[1J]. Our value of v 
obtained above is very close to this value. However, it 
should be remarked that the simple separation of vari- 
ables in terms of (p± is not perfect due to the higher- 
order terms in the compact gauge theory. Nonetheless, 
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FIG. 3: Phase diagram for N = 2. There are two phases, con- 
finement and Higgs, separated by second-order phase transi- 
tion line. There also exists a crossover line in the confinement 
phase separating dense and dilute instanton-density regions. 



our numerical studies strongly suggests that the phase 
transition for N = 2 belongs to the universality class of 
the 3D XY model. 

It is instructive to see the behavior of the instanton 
density p. We employ the definition of p in the 3D 
U(l) com pac t lattice gauge theory given by DeGrand and 
Toussaint [15j . p in FigJ5] decreases very rapidly near the 
phase transition point. This indicates that a "crossover" 
from dense to dilute instanton "phases" occurs simulta- 
neously with the phase transition. In other words, the 
observed phase transition can be interpreted as a con- 
finement (small ci)-Higgs (large c±) phase transition. 

In Figj3l we present the phase diagram for N = 2 in the 
c%-c\ plane. There exists a second-order phase transition 
line separating the confinement and the Higgs phases. 
There also exists a crossover line similar to that in the 
3D N = 1 U(l) Higgs model[ll|]- 

Let us turn to the N = 3 case. Among many pos- 
sibilities of three ci a 's, we first consider the symmetric 
case en = C12 = C13 = ci. One may expect that there 
are two (N —1 = 2) second-order transitions that may 
coincide at a certain critical point. Studying the N = 3 
case is interesting from a general viewpoint of the critical 
phenomena, i.e., whether coincidence of multiple phase 
transitions changes the order of the transition. We stud- 
ied various points in the ci — c\ plane and found that 
the order of transition changes as C2 varies. In Fig [51 we 
show E and p along C2 = 1.5 as a function of c±. Both 
quantities show hysteresis loops, which are signals of a 
first-order phase transition. In FigJSl we present C at 
c 2 = 3.0. The peak of C at around c\ ~ 0.48 develops 
as L is increased, whereas E shows no discontinuity and 
hysteresis. Therefore, we conclude that the phase transi- 
tion at (c2,ci) ~ (3.0,0.48) is second order. In FigE^a), 
we present the phase diagram of the symmetric case for 
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FIG. 4: (a) Internal energy E and (b) instanton density p for 
N = 3 at C2 — 1.5 and L — 16. Both exhibit hysteresis loops 
indicating a first-order phase transition at ci ~ 0.551. 
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FIG. 5: (a) Specific heat for N = 3 at c 2 = 3.0. (b) Close-up 
view near the peak. The peak develops as L increases. 



N = 3, where the order of transition between the confine- 
ment and Higgs phases changes from first (smaller C2) to 
second order (larger C2). In Fig^b) we present C along 
ci = 0.2, which shows a smooth nondeveloping peak, p 
decreases smoothly around this peak. These results in- 
dicate crossovers at C2 — 1.5. 

Then it becomes interesting to consider asymmetric 
cases, e.g., en 7^ C12 = c 13 . This case is closely related to 
a doped AF magnet. <j> 2 and (j> 3 correspond there to the 
CP 1 spinon field in the deep easy-plane limit, whereas 
(j>i corresponds to doped holes. This case is also relevant 
to cosmology because the order of Higgs phase transition 
in the early universe is important in the inflational cos- 
mology. Furthermore, one may naively expect that once 
a phase transition to the Higgs phase occurs at certain 
temperature T, no further phase transitions take place at 
lower T's even if the gauge field couples with other Higgs 
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FIG. 6: (a) Phase diagram for the N = 3 symmetric case. 
The phase transitions are first order in the region C2 < 2.25, 
whereas they are second order in the region C2 > 2.25. 
There exists a tricritical point at around (02, ci) ~ (2.25, 0.5). 
Crosses near C2 = 1.5 line show crossovers, (b) Specific heat 
for N — 3 at ci = 0.2. It has a system-size independent 
smooth peak at which a crossover takes place. 
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FIG. 7: (a) Specific heat of the ci = (1,2,2) model (N=S) at 
C2 = 1.0. (b,c) Close-up views of C near (b) en ~ 0.35 and 
(c) cu ~ 0.52. 



bosons. However, our investigation below will show that 
this is not the case. 

Let us consider the case C12 = C13 = 2cn, which we call 
the Ci = (1, 2, 2) model, and focus on the case C2 = 1.0. 
As shown in FigJTJa), C exhibits two peaks at Cn ~ 0.35 
and 0.52. Figs [7(b), (c) present the detailed behavior of C 
near these peaks, which show that the both peaks develop 
as L is increased. We conclude that both of these peaks 
show second-order transitions. This result is interpreted 
as the first-order phase transition in the symmetric N = 3 
model is decomposed into two second-order transitions in 
the ci = (1,2,2) model. 

Let us turn to the opposite case, C12 = C13 = 0.5cn, 
i.e., the c\ = (2, 1, 1) model at C2 = 1.0. One may expect 
that two second-order phase transitions appear as in the 
previous ci = (1, 2, 2) model. However, the result shown 
in Fig[5] indicates that there exists only one second-order 
phase transition near cn ~ 1.08. The broad and smooth 
peak near cn ~ 0.85 shows no L dependence and we 
conclude that it is a crossover. This crossover is similar 
to that in the ordinary N — 1 gauge-Higgs system as we 
shall see by the measurement of p below. 

The orders of these transitions are understood as fol- 
lows: In the ci = (1,2,2) model, as we increase cn, 
the two modes 4> xa {a = 2,3) with larger c\ a firstly be- 
come relevant and the model is effectively the symmetric 
N = 2 model. The peak in FigJTfb) is interpreted as the 
second-order peak of this model. For higher en's, the 
gauge field is negligible due to small fluctuations, and 
the effective model is the N = 1 XY model of 4> x \. It 
gives the second-order peak in Fig[7tc). Similarly, in the 
Ci = (2, 1, 1) model, cf) xl firstly becomes relevant. The ef- 
fective model is the N = 1 model, which gives the broad 
peak in Fig [5] as the crossover [llj . For higher en's, the 
effective model is the N = 2 symmetric model of <f> X 2, <f> X 3 
and U xfi , giving the sharp second-order peak in Figj8] 

In FigOS we present p of the ci = (1, 2, 2) and (2, 1, 1) 
models as a function of c\\. p of the c± = (1, 2, 2) model 
decreases very rapidly at around cn ~ 0.35, which is the 
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FIG. 8: Specific heat of the ci = (2, 1, 1) model at C2 = 1.0. 
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FIG. 9: Instanton density p at C2 = 1.0 in the (a) ci = (1, 2, 2) 
model and (b) ci = (2, 1, 1) model. 



phase transition point in lower en region. On the other 
hand, at the higher phase transition point, en ~ 0.52, p 
shows no significant changes. This observation indicates 
that the lower-cn phase transition is the confinement - 
Higgs transition, whereas the higher-cn transition is a 
charge-neutral XY-type phase transition. 

On the other hand, p of the c\ — (2, 1, 1) model de- 
creases rapidly at around en ~ 0.85, where C exhibits a 
broad peak. This indicates that the crossover from the 
dense to dilute-instanton regions occurs there just like in 
the N = 1 case 11| . No "anomalous" behavior of p is ob- 
served at the critical point en ~ 1.1, and therefore the 



phase transition is that of the neutral mode. 

We have also studied the symmetric case for N = 4, 5 
at C2 = 0. Both cases show clear signals of first-order 
transitions at d ~ 0.89(7V = 4),0.86(iV = 5). On the 
other hand, at C2 = oo, the gauge dynamics is "frozen" to 
U X fi = 1 up to gauge transformations, so there remain N- 
fold independent XY spin models, which show a second- 
order transition at c± ~ 0.46. Thus we expect a tricritical 
point for general TV > 2 at some finite C2 separating first- 
order and second-order transitions. 

Let us summarize the results. For N = 2 there is a crit- 
ical line C\{c2) of second-order transitions in the — c\ 
plane, which distinguishes the Higgs phase [c.\ > c\) and 
the confinement phase (ci < Ci). This result is consis- 
tent with Kragset et al.[lj|. For N = 3 there is a similar 
transition line, but the region < C2 < C2 C ~ 2.25 is of 
second-order transitions while the region C2 C < C2 is of 
first-order transitions. To study the mechanism of gener- 
ation of these first-order transitions, we studied the asym- 
metric cases and found two second-order transitions [in 
the ci = (1,2,2) model] or one crossover and one second- 
order phase transition [in the c\ = (2, 1, 1) model]. The 
former case implies that two simultaneous second-order 
transitions strengthen the order to generate a first-order 
transition. Chernodub et al.[16| reported a similar gen- 
eration of an enhanced first-order transition in a related 
3D Higgs model with singly and doubly charged scalar 
fields. We stress that the above change of the order is 
dynamical because (1) It depends on the value of (2) 
Related 3D models, the CP 1 *- 1 and iV-fold CP 1 gauge 
models, exhibit always second-order transitions (See the 
last reference of Ref . [2j ) ■ 
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